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Abstract 

Distance matrices are matrices whose elements are the relative dis- 
tances between points located on a certain manifold. In all cases con- 
sidered here all their eigenvalues except one are non-positive. When 
the points are uncorrelated and randomly distributed we investigate 
the average density of their eigenvalues and the structure of their 
eigenf unctions. The spectrum exhibits delocalized and strongly local- 
ized states which possess different power-law average behaviour. The 
exponents depend only on the dimensionality of the manifold. 

1 Introduction 

In a recent work about general properties of complete metric spaces 
A.M. Vershik introduced a specific type of random matrices, which he called 
distance matrices, and asked questions about their statistical properties. 

Distance matrices are defined for any metric space X with some prob- 
ability measure on it in the following way. Consider points randomly 
distributed on X according to the measure. The matrix element Mij of the 
N X N (real symmetric) distance matrix M equals the distance on X be- 
tween points i and j. In all cases considered here it is tacitly assumed that 
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there always exists a distance 1 1 . . . 1 1 on X between two points i and j which 
depends only on their relative position and we use the notation 

Mij = \ \xi- XjW, (1) 

where Xi is the c?- dimensional vector locating the point i on X, and d is the 
dimension of the base manifold. 

For any realization of the random points, (real) eigenvalues A„ and eigen- 
vectors u^"'^ of distance matrices are well defined 

N N 

Y^M^.uf = E - ^jW^f = ^nut^- (2) 

We are interested in their statistical properties. 

The first quantity to be considered is the average eigenvalue density de- 
fined as 

1 ^ 

P(A)=<^E^(A-A,)>, (3) 

or equivalently the average integrated eigenvalue density, i.e. the average 
staircase function 

N(A)=<^Ee(A-A,)>. (4) 

Here < . . . > denotes an average taken over realizations. A^N(A) (counting 
function) counts the number of eigenvalues up to the value A. 

As the elements of distance matrices are non-negative there is one large 
positive eigenvalue whose existence follows from the Perron-Frobenius the- 
orem (see e.g. V. 2, p. 49). When the metric space X is Euclidean 
or spherical all other eigenvalues have the remarkable property of being 
non-positive (the counting function associated to a distance matrix satis- 
fies N(O^) = {N — 1)/N). The proof of this for Euclidean spaces is given 
in [3] (see |1] for a general discussion of the subject). The purpose of the 
present note is the investigation of the asymptotics, in the limit of a large 
number of points N, of the average eigenvalue density at large and small 
negative eigenvalues and the properties of the corresponding eigenfunctions 
for distance matrices built from a uniform distribution of uncorrelated points 
on a base manifold. 

The plan of the paper is the following. In Section |2l one-dimensional 
spaces are considered in detail. In Section EIU we demonstrate that the case 
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of the interval is equivalent to the one-dimensional Anderson model with 
diagonal disorder. Though all states are localized, the localization length 
increases for large negative eigenvalues as discussed in Section 12 .21 When 
the localization length is much larger than the system size, the concept of 
localization becomes meaningless and a plane wave description of such states 
is adequate. This happens for large negative eigenvalues and as shown in 
Section 12.31 it leads to a power-law behaviour of the eigenvalue density. For 
small negative A, states are strongly localized and in Section IT^ it is demon- 
strated that in the one- dimensional case the eigenvalue density tends to a 
constant. The properties of the participation ratio are also discussed in this 
Section. 

When instead of the interval the circle is considered, two new phenomena 
appear. First, as demonstrated in Section El the delocalized eigenvalues 
corresponding to large negative A form quasi-doublets whose splittings are 
much smaller than the distance among them. Second, as shown in Section 
13.11 the localized eigenf unctions of the distance matrix on the circle are, 
in general, localized not in one but in two diametrically opposite regions 
(forming a kind of echo). 

In Section |3](i-dimensional spaces are investigated. First in Section l^?T] we 
introduce the continuous approximation valid for large negative eigenvalues 
and show that it leads to a |A|^'^^'^+^^/*^'^+-'^^ asymptotics of the average eigen- 
value density. In Section 14.21 it is demonstrated that if the base manifold 
has a symmetry group, large negative eigenvalues of its distance matrix form 
quasi-multiplets whose dimensions equal the dimensions of the irreducible 
representations of the group. For small |A|, the splitting of these multiplets 
becomes comparable to the distance among them and the quasi-multiplet 
structure disappears. The general condition for the applicability of the con- 
tinuous approximation is discussed in Section 14.31 where it is demonstrated 
that the quasi-multiplets are present only for the first largest negative 
eigenvalues. In Section 14.41 the behaviour of the average eigenvalue density 
for strongly localized states is investigated and it is shown that it vanishes 
as lAI"^"^. To investigate eigenfunction properties the participation ratio is 
considered in the same Section. The presence of a localization echo is also 
established for higher dimensions. 

Numerical calculations when the points are distributed uniformly on sphe- 
res and cubes of different dimensions are consistent with the results found. 
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2 One-dimensional spaces 

2.1 Distance matrices on an interval 

Let us consider N uncorrelated points xj distributed uniformly on an interval 
of length L. The distance matrix in this case is 

Mij = \xi- Xj\, (5) 

where | . . . | stands for the usual modulus. The eigenvalue equation Q reads 

N 

\xi — Xj \uj = Aui, ioT i = 1, . . . ,N. (6) 

i=i 

The eigenvalues of distance matrices are insensitive to the ordering of the 
points but the understanding of the structure of the eigenvectors depends 
heavily on it. We will rearrange the points Xj in increasing order 

< xi < X2 < . . . < < L. (7) 

Subtract Eqs. (jH} with indices i + 1 and i (assuming ((Tj)), then 

N 

A(lii+1 - Ui) = - Xj\)Uj. (8) 



As 

I II . , \ / -1 when y > Xi+i . . 



one gets 

i N 

A{ui+i-Ui) = {xi+i-Xi)[-J2'^j+ J2 ^j]- (10) 

j=l j=i+l 

After simple manipulations one proves that these equations are equivalent to 

A{Ri+i - 2R, + Ri^i) = 2(x,+i - x^)R„ (11) 

where Ri = Li — Ln/2 with Lj = Z]}=i Uj for i = 1, N and Lq = 0. 

This second order difference equation has to be completed with boundary 
conditions. The first follows from the definition of Ri 

Rn = —Ro- (12) 
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The second one can be obtained from any of Eqs. (jH)) by expressing through 
Li. Combining Eqs. (jH)) with i = 1 and i = N one gets 

A(-i?i - 2Rn + Rn-i) = 2(xi - xn)Rn- (13) 

This condition can be casted in the form of Eqs. (fTT|) by introducing the 
point xat+i = Xi. Then Eqs. (II 1|) are vahd for alH = 1, . . . , + 1 and the 
boundary conditions correspond to the anti-symmetric solutions 

Rn = —Rq, Rn+1 = —Ri- (14) 

Eqs. flll|) coincide with the one-dimensional Anderson model 

-{E~ V,)R, + = 0, for 2 = 1, . . . , AT + 1, (15) 

with diagonal disorder 

E-\/, = 2(l + |), (16) 

where ( random variables equal to the distance between 

adjacent points. When N —>■ oo and the points Xj are uncorrelated k are 
independent random variables with the Poisson distribution 

P(/) = pexp(-pO, (17) 

where 

A^ 

(18) 

is the mean density of initial points. 



2.2 Localization length 

It is well known (see e.g. that all solutions of the one-dimensional An- 
derson model ()15|1 are exponentially localized i.e. they have asymptotically 
the following decay from their maximum value, say at uq 

\Rn\ ~ e-l'^-^ol/^'-, (19) 

where koc is the dimensionless localization length. 

When |A| ^ cxD the fluctuating part of the random potential (fT^ tends to 
zero and it is convenient to use the perturbation theory developed in The 
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first terms of tlie expansion of the localization length for the one-dimensional 
Anderson model (jl5|) with a random potential eV of zero mean (< V >= 0) 
are 

^ e'^'lV^ - ^^], (20) 
hoc 8a; 

where E — 2 = e^^^x. In our case e = l/(Ap), Vi = 2(/jp — 1), and x = 2e~^^^. 
By introducing the dimensionless scaled eigenvalue 

\ = pA = -A, (21) 

one has 

This expression is valid for positive A. When A is negative the first term is 
imaginary and only the second term remains 

hoc — 4A, when A — oo. (23) 



2.3 Crystal configuration 

Though for the model (fT^ all states are formally localized, only A^ sites exist 
in our problem and, as usual for finite systems, the effect of localization can 
be ignored when the change of the wave function over the system size is small 

A^ 

- « 1. (24) 

Hoc 

For large A^ Eq. (j^Hj) indicates that states with |A| > A^/4 are practically 
delocalized and all states with smaller |A| are localized. 

For delocalized states, the fluctuating part of the potential (fTB|) is unim- 
portant. Neglecting it is equivalent to investigate the spectrum of the dis- 
tance matrix for an equally spaced points configuration 

x^ = J;^L,^oTi = l,...,N, (25) 

which we call the crystal configuration. From Eq. ()15|) it follows that for this 
configuration the Ri take an especially simple form 

Ri = aq' + bq-\ (26) 
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where q is related to the scaled eigenvalue A 



The allowed values of q (and consequently of A) are determined from the 
boundary conditions ((13)). Straightforward calculations give two equations 
for q 

+ 1 = 0, (28) 

^^^^ + 1 = + (29) 

The set of solutions of both equations (except q = —1) corresponds to eigen- 
values of the crystal distance matrix. If g is a solution then 1/g is also a 
solution and both give the same eigenvector and eigenvalue and the total 
number of solutions is A^, as it should be. 

These equations have only one real solution which gives the Perron- 
Frobenius (i.e. the largest positive) eigenvalue. All other solutions have 
the form q = exp(i0) and correspond to a negative value of A 

A = ^— . (30) 

2 sin^ 0/2 ^ ^ 

For large A^, solutions of Eqs. p8|l and ()29|1 have the form 

q = exp{u/N) (31) 

with u independent on N. The corresponding eigenvalues are 

1 2A^2 

A = o7 — ; r — ^ — TT, when A^ — oo. (32) 

2sinh2(M/2A^) ' ^ ^ 

In this limit Eq. takes the form 

cosh z = z sinh z (33) 

where z = u/2. Its unique positive real solution is z ~ 1.19968 and the 
Perron-Frobenius eigenvalue A ~ 1.6671A^^. The imaginary solutions u = 
of Eq. lead to the following asymptotics (p = 2(7m — l/(7m)) + 0(r2~^). 
Together with the solutions of Eq. (j28p (p = (2^ + l)7r the allowed values 
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of are approximately 0„ vrn, with 1 <^ n <^ N and the corresponding 
scaled eigenvalues are 

An = (34) 

The counting function N(A) of large negative eigenvalues is 

1 ^ r 

N(A) = ^ E e(A - AO - (35) 

with C = 4v^/7r. 

The N ^ oo behaviour of the crystal configuration distance matrix can 
also be obtained without the knowledge of the exact solution by noticing that 
in this limit the eigenvector components Uj can be replaced by a continuous 
function u{x) with x as in Eq. (j2Sl)- In this approximation (which we call 
the continuous approximation) the eigenvalue equation ^ takes the form 

Au(x) = — / \x — y\u(y)dy. (36) 
L Jo 

Differentiating this equation twice and taking into account that \x\" = 2S{x), 
one gets 

Au"ix) = 2—uix). (37) 
ij 

Its general solution is u{x) = ae^^ + be~^^, with A = 2N/{v'^L). Substituting 
in Eq. (jHUj) one obtains equations for a and b whose compatibility conditions 
are exactly Eqs. and (jHSJ- 

The conditions of applicability of the continuous approximation are dis- 
cussed in Section 14.31 



2.4 Strongly localized states 

Strongly localized states with small eigenvalues correspond to the configu- 
ration of points Xj in which two points, say Xi and X2 are separated by a 
distance r = |xi — a;2| much smaller than the mean distance between points 
pr <^ 1. Let Ui and U2 be the eigenvector components at points Xi and X2- 
The eigenvalue equation ^ gives 

Aui = ru2 + E l-^i ~ -^il ""j' = rui + \x2 — Xj\ uj. (38) 
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As Xi ~ X2 the sums above are approximately equal and by subtracting these 
equations one obtains that to leading order in pr 

A = -r. (39) 

Starting from this value of A it is possible to build a perturbation theory in 
higher powers of pr. 

Therefore, each time that there exists two points anomalously close to 
each other, a strongly localized state with the eigenvalue (jHIJj) is formed. The 
density of such states is equal to the probability that two uncorrected points 
are separated by a small distance r = —A. From Eq. (jl7|) it follows that at 
negative X = pA 

p(A) ~ e\ (40) 

Strictly speaking Eq. (^01) is applicable only for very small A and it indicates 
that in the one-dimensional case p(A) tends to a constant when A — > 0. 

A convenient way to distinguish between localized and delocalized states 
is to compute the participation ratio R 

When an eigenfunction is delocalized, all uj are of the same order and the par- 
ticipation ratio N increases linearly with A^, R ~ A^. For strongly localized 
states the participation ratio is independent of A^ and R ~ lioc- Therefore, 
when the number A^ of points is fixed, the participation ratio as function of 
the corresponding eigenvalue is a constant (proportional to A^) till it becomes 
equal to the localization length and ceases to depend on A^. 

On Fig.^the results for the participation ratio from numerical simulations 
are displayed and the expected behaviour is clearly seen. R is constant and 
proportional to A^ far from the origin (on the right hand side of the figure), 
and curves corresponding to different A^ coalesce to the localization length 
when approaching the origin (towards the left of the figure). 



(41) 



3 Distance matrices on the circle 

The distance between two points {i and j) on a manifold is defined as the 
length of the shortest geodesic connecting them. For later discussion, it is 
convenient to distinguish between two kinds of manifolds. We shall call a 

9 
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Figure 1: Participation ratio corresponding to a single realization for the 
unit interval smoothed over a small window of 5\ with different number 
of points. Straight line: asymptotics of the localization length, for 
comparison. 
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manifold for which there exists a single geodesic connecting two points a 
manifold of the first kind and of the second kind otherwise. The simplest 
examples of this classification are provided by the interval and the circle for 
the first and second kind respectively. 

Let us proceed to discuss the case of the circle. For a circle of radius R 
parameterized by the polar angle the distance is 

[ 27T - \{pi - (pj\, if TT < Iv^i - v^jl < 27r ^ ^ 

This equation differs from Eq. (0) and the arguments of the previous Section 
must be slightly modified, in particular for the crystal configuration consist- 
ing of N equally spaced points located at ipj = 27ij/N. In this case the 
distance matrix M takes the form 

Mij = \\ipi-^,\\ = ^f{t-j), t,j = l,...,N, (43) 
where f{k) are the integers 

f(k) = I ^' ^ ^ ^ (44) 

\ N-k, [N/2] <k<N ^ ^ 

and [x] is the integer part of x. M is therefore a circulant matrix whose 
successive rows are obtained by cyclic permutations of the first one (see e.g. 
[7j). Its eigenvectors are the Fourier harmonics 

mJ") = e^™^/^, n,j = l,...,N (45) 

with eigenvalues 

A" = f E/We'™'/""- (46) 

k=l 

As in the preceding Section it is convenient to define the scaled eigenvalues 
A as 

The sum (HBl) takes a different form for A^ even or odd. For A^ even 



1 _ (-i)n jya 
2sm (vrn/A) 4 
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and for odd 



1 - (-l)"cos(WA^) , , , 

= ■ y MH ' n = l...,iV-l, A;v = ^^. 49 

zsm (Trn/A') 4 



A„ = A^_„--(l-(-ir)-^, (50) 



The eigenvalues with n and N — n are degenerate due to the fact that both 
Uj^^ and Uj^^* = are eigenvectors of the distance matrix ()43|). When 

n/N < 1 

similar to Eq. (j34p for the crystal solution for the interval except for exact 
two-fold degeneracies for the circle. 

As for the distance matrix on the interval, the asymptotic behaviour of 
eigenvalues with n/A^ — > for the circle can be calculated by considering 

as a continuous function u{2TTj /N). In this approximation the eigenvalue 
equation reads 

NR /-stt 

Ku{lp) = / \\ip - Lp'\\u{Lp')dLp' . (51) 



27r Jo 

Taking the second derivative one gets 



NR 

Ku"{^) = —-{u{^)-u{^-T^)). (52) 

TT 



The second term appears due to the definition (|42|) of the distance on the 
circle when \ip — ip'\ > tt. The periodic solutions of this equation are e^*"''^ 
with eigenvalues given by (j50|) . Notice that for n odd the contribution to 
Eq. (j51|) from angles close to tt is the same as the contribution from small 
angles and for n even they cancel each other. 



3.1 Strongly localized states 

The behaviour of the eigenvalue density for the distance matrix on the circle is 
practically the same as on the interval (except for quasi-degenerate doublets). 
However, strongly localized eigenf unctions for the circle differ from those 
for the interval. This is in contrast to the familiar situation (e.g. for the 
Anderson model) in which strongly localized wave functions do not depend 
on the choice of boundary conditions. The origin of this difference is to be 
found in the (unusual) growth of the matrix elements with the distance. 

Let us assume that an eigenfunction is localized in a region L of the order 
of the localization length lioc with lioc -C 1 and Ui are (large) components of 
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this eigenfunction inside this region. Consider a certain point (po (measured 
from a point inside L) at a distance large in comparison with the size of the 
locahzation region. Due to locahzation the value uq of the eigenfunction at 
this point decreases exponentially |mo| ~ e~''^°/''°'='. On the other hand uq has 
to be computed from Eq. where the sum can be restricted to points lying 
in the localization region 

Auo = ^ llv^o - V^ill (53) 

Let < — v^j < TT. Then 

Xuq = (pQ^Ui-^ipiUi. (54) 

As the eigenfunction considered is a localized state, \uq\ should be much 
smaller than \ui\ for all v'o ^ hoc- But the sums on the right hand side 
include only the values of Ui inside L which are (almost) independent of 
ifQ. Therefore, in order to obey the localization property, the Wj's inside the 
localization region should satisfy 

^Mi^O, ^v^iMi^O. (55) 

Here the sign ^ means that these sums should be exponentially small 
(~ Q-M/kocy When only powers of hoc/fo are considered, the above sums 
are zero 

^ Mi = 0, ^ V9j = 0, (56) 

which can be interpreted as conditions for the vanishing of the total charge 
and the total dipole moment of charges Ui located at yjj. These are the only 
general relations to be satisfied for the interval. They do not depend on ipo 
and express the necessary conditions for the vanishing of the eigenfunction 
outside the localization region. 

However, for the circle, a new feature appears when the point y^o is close 
to a region diametrically opposite to the localization region L. In that region 
ipo = + ipQ with \iPq\ <^ 1 and, due to the definition of the distance (jl^ . 
Eq. ()53|) takes the form 

Kuo = n^Ui- ^ {ipQ - ipi)ui + ^ {i)Q-ipi)ui. (57) 
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The important difference with Eq. (|54j) is that the right hand side of this 
equation depends strongly on ip^ which determines the sphtting between 
negative and positive sums in ()57|) . No simple conditions can be imposed 
on the values of the eigenfunction inside the localization region (similar to 
Eqs. dnni)) insuring naturally that the left hand side of Eq. is small. Con- 
sequently, our assumption (required in the usual localization theory) that a 
circle eigenfunction is localized only in one small region is not correct and 
the above arguments indicate that eigenfunctions of distance matrices on the 
circle are, in general, localized not in one but in at least two diametrically 
opposite regions. Exceptions to this rule may be constituted by states local- 
ized in such a small region that the diametrically opposite one contains no 
points (i.e. one of the sums in Eq. is empty). 

On Fig 121 numerically calculated eigenfunctions of the distance matrices 
for the interval and the circle are plotted. In both figures the abscissa axis 
is the distance from the origin divided by the total length. As predicted, for 
the case of the interval (left hand side) each eigenfunction is localized in one 
small region whereas for the circle (right hand side) the eigenfunctions are 
large in two diametrically opposite regions (regions whose abscissas differ by 
a value of 1/2). 

Later we will show that the sort of 'echo' discussed here is present, in 
general, for distance matrices on manifolds of the second kind. Examples are 
given in the next Section. 

4 Higher-dimensional spaces 

In this Section we generalize the methods developed for the one- dimensional 
case to higher-dimensional spaces. 

4.1 Continuous approximation 

The asymptotics of the average eigenvalue density at large negative eigen- 
values is related to delocalized states whose contribution can be calculated 
in the continuous approximation. It is then necessary to solve the following 
equation 




(58) 
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Figure 2: Individual eigenfunctions corresponding to the n-th eigenvalue 
with N = 1000 points on the unit interval (left hand side) and on the unit 
circle (right hand side). From top to bottom: n = 200, 210, 220, 230, and 
240. 

where the points x and y belong to a dimensional base manifold X of 
volume V. 

In a small vicinity of each regular point the manifold can be considered 
as a part of the d-dimensional Euclidean space with coordinates z. In 
such vicinity eigenfunctions of Eq. can be considered as functions of z 
and we seek for semiclassical-type solutions 

u{z)r~.e''^' (59) 

with a large vector q. 

Locally Eq. (j3Hj) leads to the following expression for the eigenvalues 

H'l)-TF I W'^dz. (60) 

This formula is valid for manifolds of the first kind where two points can 
be connected by a single geodesic. For manifolds of the second kind (like 
spheres) there exist a few regions which will contribute to A(g). For clarity 
only the first case will be considered in detail. 
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Formally the integral is divergent but it can be computed from the 
convergent integral 

I(a,q)= [ e-^^^e^^dz (a > 0), (61) 

by using A(g) = —{N/V)dl{a, q)/da\a=Q- The integral can be expressed 
through the Beta function (see e.g. V. 1, 1.5.1) and the final result is 

Hq) = -^<i-iid - iy-y{qY^'^'^^\ (62) 

where ^d-i = 2tc'^^'^ /T{d/2) is the volume of the {d — l)-dimensional unit 
sphere + . . . + = 1, and r(a;) is the Gamma function. 

When the base manifold is a part of the d-dimensional Euclidean space 
and d is odd, the result (jU^ can also be obtained by successive differentiation 
of both sides of Eq. ()58p similar to what was done in Section ITTl In this case 
one also obtains exact relations between the eigenfunctions of the distance 
matrix and those of the Laplacian for d = 1, bi-Laplacian for d = 3, etc. 

For any smooth boundary conditions the density of solutions of the form 
dnni) is asymptotically the same as for the spectrum of the Laplacian (A + 
— given by 

where V is the volume of the manifold. 

From Eqs. ()62|) and the following estimate of the tail of the integrated 
density of eigenvalues of distance matrices is obtained 

N(A) ^ ^/"e(A-A(g))p(g)rfg 

- {2nYNd -^'\V) ^'^^ ' ^^^^ 

where g(A) is the inverse of the function A(g) defined in Eq. and Cd is 
a constant depending only on the dimensionality of the system. 
Introducing the dimensionless scaled eigenvalues 

\ = K^f'\ (65) 
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Figure 3: Averaged (over 50 realizations) staircase function with = 1000 
points in the unit hyper-cube of dimension d = 1,2,3,4,5 from bottom to 
top. Straight dotted hues of slope —d/{d + 1), as predicted by Eq. (jUBj) . for 
comparison. For clarity curves are shifted vertically hj d — 1 units. 

this result can be rewritten in the universal form 



where N(A) is the counting function in the variable A. 

For manifolds of the second kind like spheres the only modification of the 
above results is a slight change of the scaled eigenvalue 



where g is the number of singular regions contributing to Eq. (jHOI)- For 
spheres of arbitrary dimensions g = 2. 

On Fig. 01 results of numerical calculations of the average staircase func- 
tion for hyper-cubes of different dimensions are compared with the prediction 
fl66|) . They are in very good agreement. 

4.2 Quasi- multip lets 

The estimates of the previous Section are general but they do not take into 
account the fine structure of the eigenfunctions. 




(66) 




(67) 



17 



Let us assume that the manifold X is invariant under a certain symmetry 
group G. As the kernel ||x — ^|| of Eq. (|3H|) is the distance between two 
points on this manifold, it remains unchanged under simultaneous transfor- 
mation of X and y. Therefore from Eq. ()58|) it follows that the transformed 
eigenfunction 

u'{x) = u{G{x)) (68) 

is also a solution of this equation. From this simple remark it is clear that in 
the continuous approximation the eigenfunctions u{x) form irreducible repre- 
sentations of the symmetry group of the initial manifold similar to solutions 
of the Laplace equation. When this group has /i- dimensional irreducible rep- 
resentations the eigenvalues of distance matrices will be /i-times degenerate. 

To be specific let us consider in detail the case of the ci-dimensional sphere 
as the base manifold. The invariance group of the sphere is the d-dimensional 
rotation group. Let p = d—1. It is well known (see e.g. V. 2, 11.2) that the 
harmonic polynomials of degree p + 2 (hyper-spherical harmonics) Yirji{6, (p) 
form the basis of irreducible representations of the rotation group. Here 6 = 
9i, . . . ,9p and (p are the standard hyper-spherical angles and rh = nii, . . . , nip 
are integers obeying the inequalities 

<mp < . . . <mi <p + 2. (69) 

The dimensions h{l,p) of these representations are 

h{l,p) = {2l+p)^±^^ (70) 

(for d = 2, h{l, 1) = 2/ + 1, and ¥[^{0, ip) are the usual spherical harmonics). 

The eigenvalues corresponding to these eigenfunctions have multiplicity 
h{l,p). Their explicit form can easily be derived directly from the invariance 
of Eq. (jHH|) under rotation. Choose the z-axis along the vector x. Introducing 
the hyper-spherical coordinates in the usual way one concludes that u{y) 
equals the unique harmonic polynomial depending only on cos^ where 9 is 
the angle between vectors x and y which is proportional to the Gegenbauer 
polynomial Cf'^^(cos6') (see e.g. |SJ V. 2, 11.2). Therefore 

TV /o 

A^ = Cip— ecr{cose)smPede, (7i) 

V Jo 

with a constant Cip depending on p and /. The explicit form of is not 
instructive for our purposes. 
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Figure 4: Averaged staircase function with = 1000 points on the d- 
dimensional unit sphere (rf = 2, 3, 4, 9) showing quasi-degenerate muhiplets. 
For clarity, the counting function has been normahzed to A^. 

From properties of the Gegenbauer polynomials (see e.g. |Bj V. 2, 10.9) it 
follows that all A; with even / ^ are zero and consequently beyond the reach 
of the continuous approximation (see Section lOj) . The value corresponding 
to / = is the Perron- Frobenius eigenvalue. We therefore concentrate on odd 
values of /. 

In the continuous approximation the eigenvalues for rf- dimensional sphe- 
res are h{2k -|- 1, d — 1) times degenerate. For the one-dimensional sphere 
(i.e. the circle) h{2k + 1, 0) = 2 as has been seen in Section 01 For the 
2-sphere (the usual sphere) h{2k + 1, 1) equals 3, 7, 11 . . ., for the 3-sphere 
the first multiplicities are 4, 16, 36, . . ., for the 4-sphere they are 5, 30, 91, ... , 
etc. In general the first multiplet for the ci-sphere corresponding to / = 1 has 
multiplicity c? + 1, i.e. it is equal to the dimension of the embedded space. 
On Fig. 0] these degeneracies can be read off from the numerically calculated 
average staircase function for spheres of different dimensions. 

We also show on Fig. the structure of some individual eigenstates. We 
proceed as follows. Given an eigenvector m, choose its largest (in absolute 
value) component, say mi, and plot the components Uj as a function of the 
distance | |xi — XjH (normalizing to the maximum distance between two points 
on the manifold). By construction, the highest component corresponds to 
the zero value of the abscissa. For the 1-dimensional case, this procedure 
corresponds to the natural ordering of points on the interval or the circle. The 
quasi-multiplet structure of the eigenstates is clearly visible on the figure. For 
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Figure 5: Individual eigenfunctions u^'^^ corresponding to the n-th eigenvalue 
with = 1000 points on the unit 3-sphere (left hand side) and on the unit 
2-sphere (right hand side). From top to bottom: n = 1,5, and 21 for d = 3 
and n = 1,4, and 11 for d = 2. These values of n correspond to the lowest 
eigenvalue of each of the first three quasi-multiplets. The Perron-Frobenius 
eigenfunction {n — 1000) is at the bottom. See text for further explanation. 



comparison, the Perron-Frobenius eigenstate is also displayed. For spheres 
its components are constant to within 1/y/N fluctuations. 

4.3 Condition for applicability of the continuous ap- 
proximation 

The continuous approximation is based on the well known fact that under 
quite general conditions the sum of a large number of independent random 
variables Xj with distribution dii{x) tends to its mean value 

1 ^ r C 

E/(^.) = J fi^)dm + (72) 

where, when — cxd, C is a random variable with zero mean and variance 
independent on A^. 

This type of 'ergodic theorem' makes natural to consider, instead of the 
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true eigenvalue Eq. 



1 ^ 

its continuous approximation 

1 r 

Al^'>Un{x ) = - J^\\x-y\\ Un{y )dy (74) 

(we introduce for convenience the factor 1/(A^ — 1) with the corresponding 
redefinition A„ = A„/ (A^ — 1) because in the sum in (j73|) there are only 
— 1 non-zero terms). Here for simplicity uncorrelated points Xj uniformly 
distributed on the base manifold X (i.e. djji^x) = dx/V) are considered. 

As the kernel of the integral equation (ff^ is symmetric, its eigenvalues 
A^*^) are real and its eigenfunctions m„(x ) can be chosen real orthogonal 

^ f Unix )Um{x )dx = 6mn- (75) 

Let us first consider the case when the eigenvalue A^*^) of the continuous 
equation (ffij) is non degenerate and let us look for solutions of the true 
eigenvalue equation (ffH|) in the form 

u'f^ = Un{Xj) + J2 CmUmiXj), (76) 

where Cm are considered as small quantities. (This procedure is often used 
in perturbation theory when dealing with the Schroedinger equation.) Sub- 
stituting these expressions in Eq. (fTHj) . multiplying both sides by Un{xi), and 
summing from z = 1 to A^ one obtains 

" A*" ^ ^ ^"('^^) I ^»i('''i) ^ ' CjyiU'i^{x,i] 
1=1 \ rtiy^n 



N 



— J2 ^n{Xi) \ \Xi -XjW Un{Xj) + ^ C^Um{Xj) | . (77) 
' i,j=l \ m^n 



The mean values of the terms which are multiplied by Cm are zero and, 
consequently, they are of the order of (see Eq. ()72|l ). Therefore in 
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Eq. (|77|) these terms can be ignored to first order in 1/yN and this equation 
is reduced to 

KAn^ti) = BN{n), (78) 



where 



1 ^ 



AN{n) = —J2^niXi)Un{Xi), (79) 



i=l 

N 



N 

1 

_ n 

When N ^ oo these sums tend to their mean value plus corrections 



A,(n)-1 + ^, 5.(n) - A(^) + (81) 



The mean value of cr(ri) and S(n) is zero. When N oo their variances are 
independent on and can be computed from straightforward calculations. 

By taking the average of AnIu) and BN{n) one recovers the continuous 
approximation result. The first correction 5A„ = A„ — A^*^^ is determined by 
their fluctuating part 

5A„ = -^(E(n)-A(r)a(n)). (82) 



The mean value of (5A„ is, of course, zero and its variance is equal to 

< SAl >= -j;^v\n), (83) 

where 

v^(n) = I u'tix )dfi(x) - 1, (84) 
Jx 

provided that the Un{x ) are orthogonal. For semiclassical-type solutions 
(I59|l . where real orthogonal solutions can be chosen in the form ^/2 sin qz 
and v^cosgS*, v'^{n) = 5, independent of n. 

If the eigenvalue A^*^) of the continuous equation (ffljl is /i-fold degenerate, 
a simple extension of the previous arguments leads to an estimate of the 
splitting of the degeneracy by including the next order corrections. 
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Let ul^\x ) m = 1, . . . ,hhe the solutions of the the continuous equation 
fl7i|) corresponding to an /i-fold degenerate eigenvalue A^'^) normalized as 
follows 

/ = 6mm'. (85) 

The total splitting AA„ = EJ=i (^A^f ) is a random variable which has the 
following estimate 

^ = h.C-'-^ (86) 
with < C >= and < >= 1 and 

h h 



J^E E / [ut\x )u^:-'\x )f d^^ix) - 1. (87) 



m=lm'=l' 



More generally, if a function /(A„) is computed in the continuous approxi- 
mation, it will have a random (depending on realizations) fluctuation 5/(A„) 
which to first order in 1/ ^/N is 



We are interested mainly in the counting function. According to Eq. (jMj) it 
has a power-law asymptotics when smoothed over a small window 6A and 
dN(A)/A = z/N(A)/A with u = -{d + l)/d. Consequently 

(5N(A) ~ N(A)-^. (89) 



The condition that perturbed eigenvalues do not deviate from the ones com- 
puted in the continuous approximation (or equivalently that the splitting 
(I8H|1 is much smaller than the distance between multiplets in the continuous 
approximation) leads to the following inequality 

5N(A) « 1 (90) 

Together with ()89|) it implies that, in general, the continuous approximation 
gives correctly the eigenvalues of only the first \/N eigenstates ordered by 
increasing eigenvalues and that by increasing more and more distinct and 
isolated multiplets are present. In terms of scaled eigenvalues, one has that 
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the continuous approximation can be used for the approximation of non- 
averaged eigenvalues satisfying 

|A| > iv('^+i)/2'^. (91) 

For smaller |A| the spacing between successive eigenvalues computed in the 
continuous approximation becomes comparable to the fluctuating corrections 
and the mixing of different states is important. 

Nevertheless, one can still use the continuous approximation for the cal- 
culation of average quantities (e.g. the counting function). The main point 
is that the first correction to (jSU)) vanishes because the mean value of ( is 
zero to first order in 1/\/N. If there are no singularities (which is probably 
true for d > 2), the next order correction will be of the order of G{X)/N 
with a certain function G{X) and the conditions ()90|) on the average will be 
valid till A is of the order of 1. Therefore for averaged quantities one can 
use the continuous approximation for a finite fraction of the total number of 
eigenvalues. 



4.4 Strongly localized states 



The properties of strongly localized states can be estimated by slightly mod- 
ifying the arguments used for the one- dimensional case. 

For manifolds of the second kind (like for the circle), assume that an 
eigenfunction is localized in a small region L and Ui are (large) values of this 
eigenfunction inside this region. Choose the origin somewhere inside this 
region and consider a point Xq at large distance from it. The value uq of the 
eigenfunction at this point from Eq. (j2)) is 



A-Un 



El 



(92) 



As I |a^o| I ^ I l^il I) the right-hand side of this expression can be expanded into 
powers of Xi 

Auo = R'^Ui + '^C ■ Xi Ui + . . . , (93) 

where R = \ \xo\ \ and C = (9||xo||/(9xo- 

As I Mo I should decrease exponentially with R and each term in the above 
equation decreases as a different power of R, one concludes, similar to the 
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case of the circle, that in order to obey the locahzation property Ui should 
obey the relations 

^ Mi 0, ^ C ■ ^ 0, etc. (94) 

In regions not too close to a region diametrically opposite to the localization 
region ( = xo/R and these equations take the form of zero multipole moments 

^Mi = 0, ^XjMi,... (95) 

But close to the diametrically opposite region there are points for which the 
distance changes its form as in Eq. (j42p . For these points the derivative of | | 
will change sign and, in general, conditions (j95|) will be impossible to fulfill. 
Consequently, the eigenfunction in the diametrically opposite region cannot, 
in general, be exponentially small as required by the usual localization theory. 

The above discussion about the localization 'echo' is valid for manifolds 
of the second kind (i.e. if there exist a few geodesies connecting two points) 
e.g. for spheres of different dimensions. In Fig. [Ulwe present the structure of 
the localized eigenfunctions of the distance matrices for 2 and 3 dimensional 
spheres in the same way as it was done in Fig. |S| In all cases the eigen- 
functions have large values near the left and right ends of the interval, in 
agreement with localization on two diametrically opposite regions. Compari- 
son with Fig. illustrates also the completely different structure of extended 
and localized states (notice the difference in scale in both figures). 

To compute the density of strongly localized states let us, as above, as- 
sume that an eigenfunction is localized in a small region L and denote by Ui 
its large components inside L. (For simplicity we consider only manifolds of 
the first kind.) 

The eigenvalue equation Q gives 

Aui = ^ llxj — Xjll Uj. (96) 

Because we assume that all Ui are localized (are large) in a small region, all 
differences 11 of the order of the distance r between nearest points 

Aui^r^Uj. (97) 
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Figure 6: Individual eigenfunctions u^^^ as for Fig. except n = 800, 801, 
802, 803 for = 3 and n = 650, 651, 652, 653 for d = 2 from top to bottom. 
All states are localized. 
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As already indicated in Eq. (|94p . one of the necessary conditions of localiza- 
tion is 

Ui + Y,Uj = 0. (98) 

Therefore 

A ^ -r, (99) 

which means that the distribution of eigenvalues of localized states is approx- 
imately the same as the distribution of distances between two uncorrelated 
points randomly distributed on the manifold. For small distances the latter 
is proportional to the volume and, consequently, the staircase function at 
small negative A is 

N(A) ^ 1 - K,i-XY, (100) 

where is a constant depending on the dimensionality of the system. (In 
situations where points though uniformly distributed exhibit repulsion we 
expect that the exponent in Eq. pOOj) will be larger.) 

On Fig. [7| (same numerical data as for Fig. E)) the small A behaviour of 
the staircase function is displayed. It is clear that the estimate ()100p is in 
good agreement with the numerical simulations. 

Eq. ljlOOp together with Eq. which govern the small and large A 

behaviour of the average staircase function respectively, strongly suggest that 
N(A) is independent on the number of points when A^ —>■ oo. This is 
illustrated on Fig. |H1 . 

To get further insight on localization properties let us first remind some 
qualitative properties of the Anderson model. In one dimension all states 
are localized, irrespective of the size of disorder. In two dimensions pertur- 
bation theory cannot be applied and the localization length increases ex- 
ponentially when disorder decreases. For dimensions larger than two there 
is a localization-delocalization transition at finite strength of disorder and 
the threshold value increases with dimensionality of the system. For dis- 
tance matrices, the parameter governing the size of the disorder is the scaled 
eigenvalue A = A(A^/l^)^/'^, with increasing disorder as |A| — > 0. In Fig. Elwe 
present the numerically computed participation ratio for spheres of dimen- 
sions 2 and 3. For the 3-sphere the effect of the localization-delocalization 
transition (i.e. the sharp increase of the localization length at a finite value 
of A) is clearly visible. 
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Figure 7: Averaged fraction of eigenvalues above A with = 1000 points in 
the unit hyper-cube of dimension d = 1,2, 3, 4, 5 (from left to right). Straight 
lines of slope d as predicted by Eq. pOO|l . for comparison. 




Figure 8: Averaged staircase function for the unit 3-dimensional cube for 
different number of points. 
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Figure 9: Participation ratio corresponding to a single realization for 2- 
dimensional (left part) and 3-dimensional (right part) unit sphere smoothed 
over a small window 6X with different number N of points . The abscissa 
axis is log(— A) where A = AN^^'^ (left part) and A = AN^^^ (right part). 

5 Summary and conclusions 

Let points be distributed on a certain manifold. The matrix element 
of the distance matrix is defined as the distance between points i and j. 
Schoenberg 13] proved long time ago that for Euclidean manifolds all eigen- 
values of such matrices, except one, are non-positive. This property can be 
extended to other manifolds, in particular, for spherical manifolds treated 
here (see jl] for a general discussion). 

The study of the one-dimensional case (interval and circle) provides al- 
ready some clues of what are the ingredients governing the properties in 
higher dimensions and the properties related to the very simple crystal con- 
figuration are partially kept when 'disorder' is added. 

In general it is demonstrated that, if the points are uncorrelated and 
uniformly distributed on a base manifold, the average density of eigenvalues 
of distance matrices in the limit N ^ oo has a power-law behaviour for 
large and small negative eigenvalues, the exponent depending only on the 
dimension of the manifold (see Eqs. ()66|) and ()100|) ). 

The eigenfunctions of the distance matrices with large negative eigenval- 
ues are delocalized (for 1 and 2 dimensional manifolds the localization length 
is much larger than the system size) whereas the eigenfunctions with very 
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small negative eigenvalues are strongly localized. 

If the manifold possess a symmetry group, large negative eigenvalues form 
almost degenerate multiplets whose dimensions equal the dimensions of the 
irreducible representations of the group and the conditions for the presence 
of isolated multiplets are established. 

A distinction among base manifolds for which two points are connected 
by one or several geodesies is made. In the latter case the eigenfunctions of 
the distance matrix are, in general, localized not in one but in several regions 
(echo). For spheres of any dimension we find strongly localized states in two 
diametrically opposite regions. The understanding of the structure of the 
echo for more general spaces deserves further study. 

Strongly localized states are, by definition, mostly sensitive to local prop- 
erties of the base manifold. The existence of the echo shows, however, that 
certain global properties of the manifold are reflected on these states as well. 

What precedes is illustrated in detail by studying distance matrices of 
uncorrelated points uniformly distributed on hyper-spheres and hyper-cubes 
of different dimensions. 

The introduction of distance matrices in Ref. P was to a large extent 
motivated by the fact that they encode the metric properties of the base 
manifold. Our results show one possible way of solving the inverse problem, 
namely, the reconstruction of the initial manifold from the knowledge of the 
spectral properties of its distance matrix. We demonstrate that large negative 
eigenvalues can be approximated with 1 / a/ZV accuracy by the solutions of the 
continuous approximation Eq. ()58|) . As this integral equation is similar to 
the Laplace equation, one may conjecture that for this equation, similarly 
as for the former, the question raised by Kac "Can one hear the shape of 
a drum?" jH] can also be answered affirmatively, except probably for very 
special isospectral cases. 

When completing this work we became aware of Ref. |^ where Euclidean 
Random Matrices were introduced. Matrix elements are then a function of 
finite range of the relative distance between two points in Euclidean space 
^ij = /(I l^j ~ ^il I)- O'^^ of the main differences with respect to the present 
work is that in |^ the results strongly depend on the choice of / whereas 
here it is kept fixed f{x) = x and emphasis is put on the choice of the 
manifold. 
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